## Make PSRM Figures 2, US Map
## Andrew Reeves and Ryan T. Moore
## First: 27 January 2017
## Last: 19 January 2020

## NOTE: This file contains code used to create maps from individual
##  geolocations.  For privacy reasons, we have only included census
##  block lat / lon, which we have jittered.

## Load required libraries
library(maps)

load("figure2-USMap.Rdata")

##### Figure 2, US map: #####
png(file = "Figure-2.png", width = 960, height = 960)
map("state", fill = FALSE, plot = TRUE, col = "gray")
points(opLocMap$lonAdj, opLocMap$latAdj, 
       col = rgb(red=0, green=0, blue =0, alpha = .005),
       pch = 19, cex = .25)
dev.off()
